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(54) METHOD AND DEVICE FOR SEPARATING AND PRESUMING NOISE CONTAINED IN SPREAD 
SPECTRUM SIGNAL 



(57) When a spread spectrum signal and various 
kinds of noise arrive at a receiver in a mixed state, the 
SN ratio of the received and demodulated signal and 
the frequency band utilization efficiency of the SS sys- 
tem are improved by separating the noise and discrimi- 
natingly detecting the frequency components of the 
noise. Specifically, the received frame signal r(t) con- 
taining noise is subjected to reverse spectrum spread- 
ing by using the same spread code sequence m(t) of a 
desired station, DC components are removed from the 
reverse-spread output u(t), a respread product output is 
produced by again multiplying the output u(t) by the 

r (t)=So(t)+no(t) 



spread code sequence m(t), N analysis frequency com- 
ponents Wi which are i times the frame rate which is the 
DFT conversion analysis result of the product output, 
the noise waveform components contained in the 
received signal are presumably detected by solving 
simultaneous linear equations with N unknowns consti- 
tuted by using the set of Wi and noise frequency compo- 
nents Xj which are obtained by performing DFT 
conversion analysis only on the noise. The unknowns 
are the set of Xj. 
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Description 
TECHNICAL FIELD 

5 The present invention relates to a method and apparatus for spread spectrum communications which provide for 

enhanced SN ratio of received signals by estimating and detecting noise or similar jamming or interfering components 
mixed in spread spectrum signals during transmission. 

BACKGROUND ART 

10 

In spread spectrum communications data to be transmitted is used to modulate a sequence of spreading codes so 
that data having a spectrum of a relatively narrow band is spread over a wide frequency band; this system has excellent 
features such as small sending power per unit frequency avoidance of serious interference with other communications 
and inherent robustness against environment noise. 

is Fig. 1 0 illustrates in block form a typical configuration of a mobile communication system for spread spectrum com- 
munications over radio channels. A transmitter TX product-modulates a spreading code sequence from a sequence 
generator 1 by data b for transmission to obtain a base band transmission output St- The transmission output St thus 
obtained is used to modulate a carrier of a frequency f 0 from an oscillator 2 to obtain a send signal, which is sent over 
a radio channel to a receiver RX. Incidentally, it is customary in the art to use, as the spreading code sequence, a 

20 pseudo noise code sequence (a PN code sequence) of the same bit period length as that of the data b. The following 
description will be given of on the assumption that spread spectrum communications utilize an M-sequence now most 
widely used among the PN sequences. 

In the receiver RX the spread-spectrum modulated signal is fed via an antenna (not shown) to an amplifier 3, 
wherein it is amplified to a required level, the amplified signal is frequency-mixed with a local signal f L (=f 0 ) from a local 

25 oscillator 4 and the combined signal is extracted via a low-pass filter 5 to demodulate the spread signal of the base band 
(coherent demodulation being assumed). 

The base band spread signal and an M-sequence code, which is identical with the code used in the transmitter TX 
and generated by a sequence generator 6, are input into a multiplier 7, wherein they are correlated with each other. The 
correlated output is integrated by an integrator 8 for a one-frame period, then the integrated output of one frame is fed 

30 via a detector 9 and a coherent detector 10 to a control terminal of the sequence generator 6 to control the timing for 
generation of the M-sequence to be synchronized with the phase of the received signal. 

Fig. 11 is a schematic showing of the spectrum of a certain signal during transmission. Reference numeral 1 1 
denotes the spectrum of a spread-spectrum modulated signal and 12 the spectrum of environmental noise mixed 
therein. When demodulated (despread) by the M-sequence code in the receiver, the spread-spectrum modulated signal 

35 11, spread over a wide frequency band, becomes a narrow-band signal as shown in Fig. 12 and the environmental 
noise 12 a signal 14 spread over a wide frequency band. Hence, the spread spectrum communication system is now 
receiving attention as a communication system that essentially lessens the influence of environmental noise. 

Fig. 13 is a waveform diagram showing the correspondence between a spreading code output e M and binary infor- 
mation in a conventional direct spreading spread spectrum communication system (DS-SS). 

40 In Fig. 13, reference character b denotes binary data to be transmitted, T D the period of the data, T c its chip period 
and St its transmission waveform. Fig. 1 3 shows an example of employing the M-sequence code that has a 7-chip code 
length L as one period. In this instance, the output e^ is sent corresponding to a code "1" and its inverted version (e^f 
is sent corresponding to a code "0." Hence, the occupied band-widths of the data signal and the spread transmission 
output St are about f D =1/T D and about f C =1/T c , respectively, as given by the following equation 1 : 

45 

f c = Lf D [Eq. 1] 

Since the spread-spectrum modulated signal thus occupies a very wide frequency band, it permits reduction of 
noise power contained in frequency bands below the frequency f c to 1/(2L) and hence is robust against noise accord- 

so ingly. While in ordinary spread spectrum communication systems L»1 and an L-fold bandwidth is used, the number of 
simultaneous telephone calls N 3 is N 3 » L (approximately a fraction of the code length L) and the simultaneous transmis- 
sion capacity/Hz is (N3/L) times that of a time division multiple access system (TDMA). Consequently, the conventional 
spread spectrum communication system has a defect that its frequency band utilization efficiency of the channel is 
extremely lower than that of the TDMA system. 

55 The reason for which the number of simultaneous telephone calls Ns cannot be set so much large as compared 
with N is that it is impossible to make sufficiently small a cross correlation coefficient between an M-sequence M D 
assigned to a desired station and a different M-sequence M y assigned to another mobile station. In general, the spread 
spectrum communication (SS) system is also insufficient in suppressing dark noise, propagation noise by a delay wave 
due to mutti-path reflection during propagation and a decrease in the SN ratio due to attenuation of the signal caused 
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by fading- these factors mainly contribute to inefficient utilization of frequency by the spread spectrum communication 
system. The noise suppressing effect in the traditional DS-SS system is given by a process gain G P (Eqs. 2 and 3). 

G p = 10log 10 L fEq.2] 

L=T D /T c = F c /F D 

If the phase of each frequency component of incoming noise is completely random, the demodulated nose power 
(the output from the integrator 8 in Fig. 10) after demodulation at the receiving side is 1/(2L) with respect to the input 
noise power (the output from the LPF 5 in Fig. 10) as referred to previously. The cross correlation value C c between 
different M-sequences varies with their phases relative to each other and its mean value is given by the following equa- 
tion 4. 



c - /AT 



20 

Hence, the demodulated output based on the cross correlation value C c becomes what is called inter-station inter- 
ference noise. In the prior art, this gives rise to a problem that the SN ratio of the received demodulated signal is seri- 
ously degraded when the number of interfering stations is large. 

The present invention has been made with a view to obviating the above-mentioned defects of the conventional 

25 spread spectrum communication system. According to the present invention, when the spread spectrum signal mixed 
with the above-mentioned various noises (such as dark noise, propagation noise by the delay waves and fading and 
inter-station interference noise) is received by the receiver, the signal and the noise are separated from each other and 
noise frequency components are discriminated and detected so that the SN ratio of the received demodulated signal is 
dramatically improved over the value conventionally defined by G P or Cc thereby significantly increasing the frequency 

30 band utilization efficiency of the SS system. 

DISCLOSURE OF THE INVENTION 

To attain the above objective, the method of claim 1 for separating and detecting a noise component contained in 

35 a spread spectrum signal that is performed at a receiver of a spread spectrum communication system is characterized 
in that: a spreading code sequences are used for correspond to binary or multi-valued information,: a received frame 
signal r(t), composed of a spreading code sequence m(t) of a desired station and noise n(t), is multiplied by an analyz- 
ing code sequence g(t) to obtain the multiplied output u(t); a set of N complex frequency components [UJ (where i=0, 
1,2 N-1) attained by discrete Fourier transform analysis of the multiplied output u(t), the received frame signal with 

40 n(t)=0 is multiplied by g(t) to obtain the multiplied output y(t). a set of N complex frequency components ETJ is obtained 
by discrete transfer analysis of y(t); the set of [UJ of the complex frequency components is divided into the in-phase 
component [U|J and the quadrature phase component [U QI ] to the set prj; the received frame signal r(t) is further ana- 
lyzed to obtain a set [RJ of complex frequency components; the amplitude of a component R* in the set of the complex 
frequency component [RJ, which has an excessively large or small amplitude as compared with a preset correcting 

45 spectrum curve, is corrected within a range in which a value, obtained by selectively adding together the powers of the 
set [UJ and the component [UqJ (where i=0,2,... f N-1) with respect to i, decreases or remains unchanged; and the 
amplitude of a component R k , of a second excessively large or small amplitude is corrected within a range in which the 
power sum decreases or remains unchanged; these operations are repeated to reduce the power sum toward zero, 
then noise components \ and X' k are estimated from corrected values of the components Rj< and R k , and the esti- 

so mated noise components are removed from the received frame signal r(t). 

According to claim 2, the method of claim 1 for separating and detecting a noise component contained in a spread 
spectrum signal in a spread spectrum communication system, is characterized in that: of the frequency components of 
the received frame signal only the component projected to the phase of each frequency component [CJ obtained by a 
discrete Fourier transform analysis of the spreading code sequence m(t) of the desired station is extracted; and the time 

55 signal obtained by the inverse discrete Fourier transform analysis is used as the received frame signal r(t). 

According to claim 3, the method of claim 1 for separating and detecting a noise component contained in a spread 
spectrum signal in a spread spectrum communication system, is characterized in that a square wave spreading code 
sequence m(t) of the desired station or a sequence m L (t) obtained by band limiting it is used as an analyzing code 
sequence g(t). 
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According to claim 4, the method of claim 1 for separating and detecting a noise component contained in a spread 
spectrum signal in a spread spectrum communication system, is characterized in that the analyzing code sequence g(t) 
is selected so that the multiplied output of the analyzing code sequence g(t) and the spreading code sequence m(t) of 
the desired station becomes a square wave that has a period kTc (where k«..., 1/3.1/2,1,2,3,...). 

5 According to claim 5, the method of claim 1 for separating and detecting a noise component contained in a spread 

spectrum signal in a spread spectrum communication system, is characterized in that the analyzing code sequence g(t) 
is selected so that the multiplied output of the analysis code sequence g(t) and a sequence m L (t) obtained by band lim- 
iting a square wave spreading code sequence m(t) of the desired station is composed of a sine wave of a period kTc 
and a limited number of its high-frequency waves. 

io According to claim 6. the method of claim 1 for separating and detecting a noise component contained in a spread 
spectrum signal, is characterized by the use of a received signal containing noise and a time wave form obtained by 
limiting a spreading code sequence for use at the receiving side to a passband frequency band lower than a predeter- 
mined upper limit frequency. 

According to claim 7, the method of each of claims 1 through 5 for separating and detecting a noise component 

15 contained in a spread spectrum signal in a spread spectrum communication system, is characterized in that the multi- 
plication of the received frame signal or multiplied output u(t) and the receiving side spreading code sequence m(t) in 
the time domain is performed by a convolution in the frequency domain. 

According to claim 8, the method of each of claims 1 through 5 for separating and detecting a noise component 
contained in a spread spectrum signal in a spread spectrum communication system, is characterized in that: a trans- 

20 mitting signal is a time waveform obtained by the convolution of an impulse train produced by q-point sampling per chip 
time width of a square wave M -sequence and a sampling function waveform having its pass band limited to q/2; and the 
receiver uses, as the despreading signal the same waveform as the time waveform. 

An apparatus of claim 9 for a receiver demodulating system in a spread spectrum communication system in which 
a spreading code sequence is made to correspond to binary or multi-valued information, is characterized by: means for 

25 despreading a received frame signal r(t) containing noise and the same spreading code sequence m(t) as that of a 
desired station; means for multiplying again the despread AC output, obtained by removing the DC component from the 
despread output u(t), with the spreading code sequence m(t) to obtain the re-spread output w(t); means for obtaining N 
respread frequency components Wi with i-times (where i=0,1,2,...,N-1) the frame rate that is the output obtained by a 
discrete Fourier transform analysis of the respread output w(t); means for solving simultaneous linear equation with N 

30 unknown variables consisting of a set of j-th (where j=0,1 ,2,...,N-1) noise frequency components Xj obtained by a dis- 
crete Fourier transform analysis of the received frame signal with the noise alone, and the coefficient matrix which 
relates the noise frequency components Xj and the respread frequency components Wi obtained in advance, and the 
set of Wj; and means for detecting the signal component in the frame obtained by removing the estimated noise wave 
form, inverse discrete fourier transform of the set X jp from the received frame signal r(t). 

35 According to the present invention described above in brief, the transmitter produces and sends a frame signal 
which is made of a spreading code whose polarity corresponds to the information to be sent, and the receiver first ana- 
lyzes the received frame signal to estimate the noise-frequency components contained in the received frame signal, 
then removes the estimated noise from the original received signal, and performs ordinary desrpeading and demodu- 
lation processing-this enhances the SN ratio of the received and demodulated signal. 

40 Based on the fact that the despread output is obtained by despreading (multiplying) the received frame signal with 
the spreading code corresponding to the desired station and the despread AC output is obtained by removing the DC 
component from the despread output, the despread AC output is composed of an AC noise component and a noise- 
dependent DC component resulting from the despreading process, it is possible to separate the AC noise component 
and the noise-dependent DC component. 

45 Since the frequency component, obtained by a DFT analysis of the re-spread output is composed of the DC and 
AC components of the incoming noise and the noise-dependent DC component, the received noise component can be 
separated and estimated by solving simultaneous linear equation with N unknown variables, where N corresponds to 
the spreading code lengths. 

As one of the received signal, let us consider the DC component of a pseudo noise and its single-frequency com- 
50 ponent with a known power whose phase matches that of a frequency component contained in the despreading code. 
By observing the re-spread output produced by the method stated above, the coefficients of the simultaneous equa- 
tions are obtained, which can be used to solve the equations. 

By only using noise components in-phase with the despreading code among all the frequency components of the 
received noise, the processing can be simplified. 
55 The present invention is also effective against interference noise which is composed of spreading codes of the 
other stations. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of the demodulator to be used at a receiver in a spread spectrum signal communication 
system according to the present invention; Fig. 2 is a diagram of the single side spectra of an M-sequence; Fig. 3 is a 

5 block diagram illustrating an example of a pre-f itter for use in the spread spectrum signal demodulating system accord- 
ing to the present invention; Fig. 4 is a block diagram of a noise discriminating and separating circuit for use in the 
spread spectrum signal demodulating system according to the present invention; Fig. 5 is a diagram showing an exam- 
ple spectrum of a respread output for a single frequency input; Figs. 6(a) and (b) diagrams showing spectra of a 
received frame signal and estimated noise; Figs. 7(a). (b), (c) and (d) are diagrams for explaining steps of correcting the 

10 estimated noise spectrum; Figs. 8(a) and (b) are diagrams showing multiplied-time waveforms; Figs. 9(a) and (b) are 
diagrams showing a multiplied-time waveform and its spectrum; Rg. 10 is a block diagram illustrating a general config- 
uration of a spread spectrum communications system for mobile communication; Fig. 1 1 is a diagram schematically 
showing the spectrum of a signal during transmission; Fig. 1 2 is a diagram showing the despread spectrum of the signal 
and noise observed at a receiver of a spread spectrum communication system; and Fig. 13 is a diagram showing the 

15 correspondence between a spreading code and binary information in a conventional direct spreading type spread spec- 
trum communication system. 

BEST MODE FOR CARRYING OUT THE INVENTION 

20 The present invention will hereinafter be described in detail based on its embodiments shown in the accompanying 
drawings. 

Fig. 1 illustrates in block form a demodulator according to the present invention. In Fig. 1 a received input signal 
ro(t) into a pre-f ilter (FIL) is a time waveform of a frame period T D resulting from the addition of noise n 0 (t) to a transmit- 
ted signal So(t) and it can be expressed by the following equation 5. The pre-filter FIL removes from the noise compo- 
25 nent n 0 (t) added during transmission a quadrature component no(t) orthogonal to the transmitted signal s 0 (t), leaving 
only a component n,(t) in-phase with the signal s 0 (t). That is. the input noise n 0 (t) is given by the following equation 5. 

R 0 W = s 0 (t)+n 0 (t) (0*t£T D ) [Eq. 5] 

30 n 0 (t) = n,(t)+n Q (t) [Eq. 6] 

As is evident from Eq. 6, the input noise n 0 (t) is composed of two components but only the component n|(t) is left. 
The orthogonality relation between the components n, and n Q means that phase angles of respective frequency com- 
ponents, obtained by a discrete Fourier transform (DFT) analysis of the components nj and n Q . are orthogonal to each 
35 other. 

Further, the pre-filter FIL has a function of removing a high-frequency component in the input signal r 0 (t) (usually, 
a frequency component above the chip rate f c ). Hence, the filter output can be expressed by the following equation 7. 

r(t)=s 0 (t)+n,(t) [Eq.7] 
40 m s(t)+n(t) (O^t^TD) 



where s(t) and n(t) are components that are obtained by band-limiting the signal So(t) and the noise n 0 (t) and extracting 
from the noise component n 0 (t) the component n,(t) in-phase with the M-sequence (which will hereinafter be assumed 
45 to indicate the sequence of a desired station). 

The filter output r(t) is applied to a despreading noise detector (DEND) to calculate therein an in-band noise spec- 
trum estimated value X(f)(<f c ). By inverting the estimated value X(f) in an inverse DFT circuit (IDFT), its time waveform 
n*(t) is obtained. The estimated noise time waveform is fed to a subtracter wherein it is subtracted from the filter output 
r(t) to obtain an output s* given by the following equations 8 and 9. 

50 

s*(t) = r(t)-n*(t) [Eq.8] 
= S(t)+e(t) 

s(t) = n(t)-n*(t) [Eq. 9] 

55 

where e(t) is the time waveform of a noise estimation error. If a substantially correct noise estimation can be made so 
that e(f)«n(0 . the SN ratio of the output s* remarkably increases as compared with the SN ratio of the filter output r(t). 
The output s* is fed to a demodulator of an ordinary DS-SS system (corresponding to 7, 8 and 9 in Fig. 10), wherein 



5 



BNSDOCID: <EP 079B871A1J_> 



EP 0 798 871 A1 



10 



15 



SO 



25 



30 



received information b f corresponding to the transmitted one b is demodulated and detected. The demodulating circuit 
comprises a matched filter and a detector. Reference character s(t) denotes a reference signal that is applied to the 
matched filter and e, the output of the integrator which is a part of the matched filter. 

Now. a description will be given for cases where as the spreading code sequence the M-sequence shown in Fig. 
13 is used. In the first place, a time waveform mss{t) of a non-band-limited square wave M-sequence will be considered. 

The M-sequence time waveform mco(t) can be expressed by the following equation 10, letting the number of sam- 
ples per chip be represented by q. 

Incidentally, » means a suffix «>, but in this specification it is expressed in the form of » at some places due to lim- 
itations on the kinds of letters usable as suffixes (In the drawings it is represented as the suffix »). 

qL-1 

mco(t)= £ c [MV P TC (T-kT c ) [Eq. 10] 

k=0 

where T c is a chip period, Q a Gaussian symbol, q the number of samples/T c , c^' a voltage amplitude at a k-th sample 
point which is +1 or -1 V, and P-rc(t-kT c ) an isolated squared pulse waveform whose voltage is 1 V from t=kT~ to 
(k+1)T c . 

By a DFT analysis of the wave form m^f), a both side frequency component C - corresponding to an i-th frequency 
(i*D. f d~ t d ) is 9»ven by the following equations 1 1 through 13. 

1 qL-1 _j2nik 

C ' = 51 E c /*/q/ e ^ ■°/+*# [Eq. 11] 

C,-^of+pf>0 (h*0) [Eq. 12] 



[Eq. 13] 



where a b ft are real and imaginary parts Cj the amplitude value of the i-th one-side frequency component and 0j its 
phase angle.Fig. 2 shows parts of C jp a { and p j( calculated from the waveform ^ of Fig. 13 with q set at 10. (The actual 
both side spectra exist to i=Lq-1=69 .) 

On the other hand, a similar analysis of the noise n(t) gives an i-th both-side frequency component X*j by the follow- 
40 ing equations 14 through 17. 



1 qL-1 .j2nik 



[Eq. 14] 



k~0 



45 



e , i f =tan * r ^±e, 
a/ 



[Eq. 15] 



50 



X, = (sign e',)^*/ +P/ (*0) 



[Eq. 16] 



55 



X, = {signQ' jjal+fy (/=0) 



qL-1 
U0 



[Eq. 17] 
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where Xj is the amplitude of an i-th one-side frequency component, which takes a positive or negative value. (This is 
defined on the basis of the phase of the both-side frequency component CV) 

As indicated by Eq. 7, the noise component to be noted here is the component nj(t) that matches the phase angle 
Oj and the noise component n Q (t) orthogonal thereto does not appear at the output of an ordinary demodulator, and 
hence it can be neglected. The frequency component Xj provided from the pre-filter FIL in Fig. 1 is in-phase with the 
frequency component CV Therefore, letting the phase of the frequency component X*j be represented by 6'j, the follow- 
ing equation 18 holds true for Eqs. 13 and 17. 

ev = ±0 [ [Eq. 18] 



10 



Fig. 3 is a block diagram illustrating an example of the pre-filter FIL described previously with respect to Fig. 1, 
which is shown to employ a DSP (Digital Signal Processing) circuit configuration. An input frame signal r 0 (t) is converted 
by an analog-to-digital converter A/D to digital form, thereafter being subjected to discrete Fourier transform processing 
by a DFT{f c -N Q (f)}. 

is As a result, the frequency components higher than the frequency f c are removed, the noise component N Q (t) 
whose phase is orthogonal to that of each frequency component of the M-sequence is also removed. This DFT output 
R(if D )=R j is subjected to inverse discrete Fourier transform processing by an IDFT means to obtain the time waveform 
r(t) given by Eq. 7. The output R s is composed of a signal component Sj and a noise component Xj as indicated by the 
following equation 19. 



20 



40 



50 



R^Sj+Xj [Eq.19] 



Fig. 4 illustrates in block form the noise discriminator/separator (DEND) that forms the principal part of the present 
invention. The binary information for each frame set at the sending side is a constant that is 1 or -1 (0^t<T D ), but since 
25 the amplitude of the signal component at the receiving side has undergone variations during transmission, the signal 
component s(t) in the input r(t) can be expressed by the following equation using a real value b. 

s 0 (t) = bmao{t) (0<rr<TD) [Eq. 20] 

30 The actual signal component s(t) is band limited in advance by a filter at the sending side and is further band limited by 
the filter (FIL) at the receiving side, but here the ideal waveform given by Eq. 20 is considered as the input and let it be 
assumed that only noise is band limited. In such an instance, the input r(t) is multiplied by moo(t) to obtain an output u(t), 
the DC component included in this output u(t) is removed by a capacitor C, and the capacitor output is multiplied again 
by m»(t) to obtain an output w(t). Since moo(f)x m<x>( t)= 1 , the output u(t) is given by the following equation 21. 

35 

u(t) = [bm<x>{t)+n(t)]m<x>(t) [Eq. 21] 

= b+n\t) moo(t) 

m b+U D +U H {t) 

where b is a DC component corresponding to the signal s(t) and u D and u H (t) are composed of a DC deviation U 0 and 
a high-frequency component, both are produced by the noise n(t) and are given by the following equations 22 through 
25. 

45 r L-1 -, 

u D =U 0 = FA X (CVX',) [Eq.22] 

L/=(l-T) J 



L-1 

U(p)= £ {C'^X't) [Eq.23] 

-i = qL-i [Eq.24] 
p-i = qL+p-i 
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[Eq. 25] 

hwi J 



to where R is a real part and F* an IDFT symbol. 

After the removal of the DC component (b+u D ) from the output u(t), the resultant component is fed to the next stage, 
where it is multiplied again by mco(t) to obtain a re-spread output w(t). 

As will be seen from the Fig. 4 configuration, the output w(t) is obtained by multiplying n(t) with mrotf) twice and then 
subtracting u D moo(t) from the multiplied output. This relationship can be given by the following equation 26 by removing 
1 s the DC component from the signal component of Eq. 20 and then multiplying K by mco(t). 

w(t) = {u(t)-b-u D }m ^(t) [Eq. 26] 

= {n(t)m oo (t)-u D }m 0O (t) 
= n(t)m <c (t)m o0 (t)-u D m co (t) 
20 =n(t)-u D m 00 (t) 



The term of the Dc component b generated by the signal component bmm(t) in the input frame signal is removed 
because H is not an object to be detected in this instance. Since the influence of the DC component b can thus be 
25 excluded with accuracy, noise can be detected by the method described below. The frequency component in Eq. 26 is 
expressed by the following equation 27 using a DFT symbol F. 

W=F[w(t)]=F[n(t)]-u D F[ moo (t)] 

-ZPWrj [Eq271 



30 



35 



qL-1 v qL-1 L-1 
>0 j-IL-L+1 J j=0 i*.(L-1) 



Now, using the component u D given by Eq. 22 and taking into account that the components X'j and C'j are in-phase 
with each ether, W(f) can be expressed as the sum of the amplitudes Wj of i-th frequency components. The calculation 
40 of the amplitude Wj gives the following equations 28 and 29. 

L-1 

W= £W/ [Eq 28] 



45 

w 



50 



= C,C y U=0)~ 



where aq is a DC-dependent deviation coefficient which is a coefficient representing the ratio of the components Wj to 
55 the noise component Xj. That is, it corresponds to the i-th frequency component of the M-sequence which is contained 
in the output w(t), as much as the noise-dependent DC component that the noise component X, generates. 

The noise component Xj contained in the input frame signal r(t) is band limited, to frequency regions less than f c 
and hence, it results in; <,(L-1). Hence, the application of Eq. 28 to L noise spectra X, Q= 0, 1 , 2 .... L-1), inclusive of the 
DC components, gives the following simultaneous linear equations with L unknown variables in Xj: 
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[Eq. 30] 
-a^ -a^ (l-a* 2 ) 



Eq. 30 can be abbreviated as follows: 



r x. ■ 




w 0 - 


x, 










w, 


• 




• 


• 




• 









[A][X]=[W] 



{Eq- 31] 



The above coefficient matrix can be pre-calculated by using the M -sequence spectrum (Eq. 11) used. The value W, 
can be obtained as by multiplying the input r(t) by mooft), removing its DC component, multiplying the remainder by Moo 
(t) and subjecting the multiplied output to a DFT analysis. Accordingly, Eq. 31 can be solved, with the result that the 
value Xj is obtained. 

In this process, e in Eq. 8 can be set at zero theoretically. In the case where a pattern (magnitude and polarity) 
regarding j of the input noise Xj closely resembles the pattern of the value Ci, however, U H (t)=0 in Eq. 21 , so that the 
component w(t) becomes zero in Eq. 26 and hence cannot be detected. 

Next, a description will be given of means for calculating each term of [A] in Eq. 31 . 

Now consider the case of applying a single frequency component (Eq. 32) as the input signal in Fig. 4. 



nj(t) = Xj cos(27crf D t+6j) (0^T D ) 



[Eq. 32] 



where 0j indicates the phase of the j-th frequency component of the M-sequence. Letting W*| represent the i-th fre- 
quency component of the value W(f) that is obtained by simulation or measurement when X ( =1 , Eq. 33 can be 
obtained. 



W*i = -a y 
o 1-a „ 



[Eq. 33] 



0=0 



40 



45 



50 



55 



In Fig. 5 there is shown an example of the simulated output of W*j when j=1 . 

The following method can also be used to determine the value ay. In practice, Eq. 20 cannot be assumed and the 
signal components are also band limited usually as shown in Eq. 7. The input signal in such an instance is expressed 
by Eq. 34. 



s(t) = bm L (t) (o<sfsf c ) 
r L (t) = s L (t)+n(t) 



[Eq. 34] 
[Eq. 35] 



where m L (t) and m H (t) are waveforms composed solely of those of the frequency components of the waveform mooft) 
which are below and above f c , respectively. 

The spectra of signals r L (t) and s L (t) match the value in Eq. 19 but are limited to i^(L-1). In this case, the influence 
of the signal dependent component b cannot completely be excluded from w(t) by the means of the multiplication and 
DC cutoff. If the above-mentioned scheme is used, the ranks of the coefficient matrix of Eq. 30 drops from L to L-1 . This 
means that one of the simultaneous linear equations with L unknown variables is expressed as a linear combination of 
the other (L-1) equations. Hence, no solution can be obtained unless any one of the values in a set of L unknown vari- 
ables , [Xj] (j as o, 1, 2, N-1), is given in advance. Alternatively, it can be solved simply by assuming the value b. 

That is, one value X|< in the set is assumed as X*k or the value b as b* and then an estimated value [X*j] of the 
set [Xj] is calculated by the above-mentioned method. Next, it is necessary to evaluate the validity or appropriateness 
of the assumed value Xk or b* on the basis of the set [X*j]. 
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Now, letting it be assumed, for simplicity, that k=L-1 is assumed, the following equation is reduced that the ranks 
of which is decreased by one. 

[Eq. 36] 





(l-*oo) -So i -«« — 


• 


x% 




' w 0 ' 


10 


-a l0 (l-a u ) -a l2 — 








w, 








x% 








• • • • 


• 


• 




• 


15 


m m • • 


• 


• 




• 




~ *L~1J* ~ a L-Xt "~ *v-v 


(l-*L-2X-2). 









20 With the use of Eq. 36, a pre-assumed value X* L1 and the set [X*j] Q = 0, 1 , 2..., L-2) are determined. 

Xj< is assumed as X* k , using a set [RJ of one side spectra of the received input rL(t) of Eq. (35) obtained by DFT 
processing. In such an instance, the following equation 37 holds for the k-th frequency component. 



25 



30 



35 



40 



b* = (R k -X* k )/C k 



[Eq. 37] 



Thus, a set of solutions [X*J calculated by Eq. 37 assuming X* k and a similar set of solutions [X*J calculated assum- 
ing b* that bear the relationship of Eq. 37 match each other. 

The following equation shows the relationship between the solution [X*J of Eq. 36, obtained by assuming b* or X* k , 
and the received signal spectrum [RJ of Eq. 19. (In the following description, i<L unless otherwise indicated.) 



R^bCj+Xj 
X* j = X | AbC j 
Ab « b*-b 



[Eq. 38] 
[Eq. 39] 
[Eq. 40] 



In Figs. 6(a) and (b) there are shown the set [RJ, the correct value [Xj] and the estimated value [X*j] in the case where 
a noise component is present at the j-th frequency alone. 

On the other hand, the product outputs by multiplying r L (t) by mco(t) and n\(t) by moo(t) are given as follows: 



u L (t) = {m L (t)+n(t)}moo(t) 
Y(t) = m L (t)xmco(t) 



[Eq. 41] 
[Eq. 42] 



45 Sets of one-side spectra obtained by DFT analyses of the product outputs u L (t) and y(t)are represented by [Up] and[Tp], 
respectively. The value [Up] in this case differs by [Tp] from the value of the one-side spectrum calculated by Eq. 23. 

Here, letting the phase of [Tp] be represented by yp, U p is split into a component U|p in-phase with y p and a com- 
ponent U Qp orthogonal to y p . That is, if y p is used as a reference phase, U p is given as follows: 



50 



U p = U lp+ jU Qp 



[Eq- 43] 



The component [UopJ has nothing to do with the signal component bm L (t) contained in the received signal r(t). Accord- 
ingly, even rf Cj varies in Fig. 6. the estimated noise spectrum X*j changes from its true vaiue Xj by AbCj but the set [UopJ 
is free from the influence of Ab. Then, b* is set at 0 and Rj in the following equation is regarded as a primary estimated 
55 value of noise. 



X*j = Rj [Eq. 44] 

Now, a description will be given on the assumption that the value Rj shown in Fig. 7 is provided as an input. In Fig. 
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7, b m Cj is a curve proportional to an m-sequence spectrum passing through the maximum value R 4 among those Rj<0 
and b n Cj (where b n <bm) is a corrected (spectrum) curve that passes through the second largest value R 2 . Now, let it be 
assume that spectra above b n Cj are all noise components, 



5 R 4 > b nC 4 [Eq.45] 

Then, assume that the power P Q (Eq. 46) of the quadrature component diminishes or remains unchanged as the result 
of Ai reduction of the amplitude R 4 . In such a case, a slight reduction of b n and another A 2 reduction of the amplitude 
R 4 are repeated, and the power P Q is evaluated each time-this processing is continued until the power P Q reached the 
io minimum value or changes to increase. 

In this process the total amount of reduction AR 4 given by the following equation becomes an estimated value X* 4 . 



75 



2(L-t) 

p=1 



Pq= Z U op 2 



X* 4 e AR 4 = A 1 +A 2 +» 



[Eq46] 
[Eq. 47] 



20 In general, when b h is further reduced, R 2 becomes an excessively large spectrum component, then the exces- 
sively large component of the value R 2 (the portion above the broken line bj is removed, and the value of the power 
P Q is evaluated. If the power P Q increases, b h +R 2 is returned to the original value. That is, the same processing is per- 
formed immediately before the power P Q increases. (In this instance, when b n further decreases, it is assumed that the 
power P Q has increased and the value R 2 is not reduced.) Moreover, to cope with the case of Xj<0, b* h (a negative 

25 value) in Fig. 7 is assumed and the same processing as mentioned above is carried out. Fig. 7(b) shows the results of 
the reduction of the amplitudes R-j and R 4 , using b h and b* h in Fig. 7(a). Figs. 7(b), (c) and (d) show an example in which 
it is assumed that the power P Q continues to decrease (or remain unchanged) as the curve of the negative side is fur- 
ther diminished from b h to b" h and hence to b"' h as shown. In this example, in the process from Fig. 7(a) to (d), the value 
Ri once becomes zero and then increases up to R-| ,w of the opposite polarity This reveals that the noise estimated value 

30 of the first frequency represented by X* ., =R 1 -R ,M ., , is a negative large value. In Fig. 7, b^C\, b 2 Cj, ... are fluctuations 
of the corrected curve.Ri\ R-i", Ri w are values obtained by subtracting the estimated noise from or adding it to the first 
or initial value R 1 , and if the estimation accuracy is high, only the signal component tends to be gradually left remaining. 
By changing the amplitude b h little by little and performing the above processing for Rj, R'j, ... which are larger or smaller 
than b h Cj, the estimated value [X*j] is obtained. 

35 A description will be given of enhancement of the correction accuracy, with the power of the in-phase component 
represented by P t of the following equation. 

[Eq. 48] 

40 2<t-U 



45 

Since a change in the noise component Xj and an increase or decrease in the power P| bear a functional relation 
also regarding j, the estimation accuracy is increased by checking variations in the power Pj while referring to the func- 
tional relation and by adjusting the amount of AR, to be decreased or increased accordingly. A value including the DC 
component may also be used as P ( . By correcting the value Rj while checking the value P ( , it is possible to perform the 

so necessary processing without reducing the component b* of the received signal more than required. 

Furthermore, a detailed estimation can be made by a method that takes into account a pattern concerning i in U (i 
which is a component of P| given by Eq. 48. That is. if a pattern of the component U H , which corresponds to the signal 
component y(t)=m L (t)moo(t) in Eq. 41, is set as a signal pattern, the component Ujj can be split into a component [U Mi ] 
that correlates with this pattern and a component [U| Qi ] that intersects it orthogonally. Accordingly, letting that power 

55 component of the power P t corresponding to [U| Qi ] be represented by P| Q , the sum P' Q of the powers P Q and Piq is used 
as a measure for evaluation and the value Rj can be adjusted by changing the wavy line of b h in the direction in which 
to decrease the evaluation measure, (in this instance, the correction can be made by using b h on the plus side and b' h 
on the minus side alternately with each other and using the value on the side made effective by a decrease in any one 
of Pq*.) Generally speaking, this method means that the sum of values obtained by selectively multiplying the compo- 
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nerrts U\ t and U Qj by weighting coefficients is used as the evaluation measure. 

The waveform by which the received input signal r(t) is multiplied is not limited specifically to the waveform moo(t) 
but may also be a different waveform g(t) and, also in this case, the same correction processing as described above can 
be performed as exemplified in Fig. 8. As shown in Fig. 8(a), the input waveform m(t) is multiplied by g(t) to obtain a 
5 square wave y 0 (t) that has a 2T C period and a 50% duty factor. Then, a waveform y* 0 (t-T) inverted from the waveform 
y 0 (t) is added to the latter in the subsequent T sec to generate a waveform y(t) that is given by the following equation. 

Y<t) - YoW{H(t)-H(t-T)] +y * o(t-T)[H(t-T).H(t.2T)] [Eq. 49] 

10 where H(t) is a step function. y(t) includes a 7-cycle square wave that has a 2T period. Hence, by making a DFT analysis 
over the 2T period with the fundamental wave set at f * D =(2T) 1 =0.5f D , the spectrum in the case of the input being m(t) 
becomes r h indicated by the solid line in Fig. 8(b). Here, 

h = (2m-1)L (h=1,2,...) [Eq. 50] 

15 

By generating the waveform y(t) from that y 0 (t), the spectrum r h can be simplified. 

The output obtained by multiplying the received input r(t) by g(t) as is the case with Eq. (20) is defined by the fol- 
lowing equation: 

20 u(t) = [bm(t)+n(t)]g(t) {Eq. 51] 

= by(t)+v(t) 

Let Uj represent the i-th spectrum that is obtained by DFT analysis for the 2T-frame output u(t) and Vj the i-th spec- 
trum for the component v(t) corresponding to the noise. Since r h * 0 for the frequency component i=h , it is possfole to 
25 define components V H and V Qi that are in-phase with and orthogonal to the spectrum Irrespectively, with respect to 
the noise component Vj. Therefore, the spectrum Uj is given by the following equation: 

[Eq. 52] 

30 

U, = T-Vu + V Qi (' = A)j 

u, = r, {i+h)\ 

35 

In such an instance, P x can be defined by the following equation as the multiplied output power corresponding to the 
noise alone. 



Px-Z^i+E V 2 ^ [Eq.53] 

i=0 m=1 



The range of integration of the above equation is, in practice, a period in which main spectra exist In Figs. 6 and 7 
45 the i-th spectrum of the input signal has been modified using the fact how the power, P Q , of the components of the mul- 
tiplied output u(t) which is orthogonal to n(t) decreases or is unnecessary as the measure for evaluation. But in this 
example the modification is made using the multiplied output power P x or P Q *=P X +P 1Q as mentioned previously 
instead of using the power P Q . 

In general, the cycle period of the waveform y(t) can be set at T p =kT c (where k=..., 1/3,1/2,1,2,3, ...). Fig. 8 shows 
so- the case where k=2. Figs. 9(a) and (b) show the cases where k=4 and k=1, respectively. For k=4, a 3.5-cycle square 
wave with period of 4-Tc can be made by adding the square wave y 0 (t) so modified as making the time axis inverted to 
the second time frame (T-2T). 

If a 4-period waveform y(t) is produced by further adding a polarity-inverted version of the waveform of the two peri- 
ods, the square wave contains seven cycles. By a DFT analysis of the 4-period wave form, a simple spectrum similar 
55 to the previously-mentioned r h can be obtained. When k=l , the original period of the square wave y(t) contains an inte- 
gral number of cycles, and hence the DFT analysis can be conducted with the original period of the square wave y(t). 
Also when k£1 , the analysis can similarly be carried out. In the case of a prime number such as k=3, the spectrum com- 
ponent of the square wave y(t) can similarly be simplified by adding after the original period a plurality of frames gener- 
ated based on the original period as in the above. 
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Where the signal m(t) is band-limited, such a deformed waveform as indicated by the broken lines in Fig. 9(a) is 
resulted, but in this case, the waveform g(t) is also deformed as indicated by the broken lines. In such an instance, the 
waveform y(t) resembles a sine wave and high-order components of the value r h decrease and its spectrum becomes 
simpler, allowing more ease in the evaluation by the multiplied output power P x . 

5 Incidentally, letting a*y represent a component that is generated in the spectrum Uj of Eq. 51 based on the noise 

component Xj of unit power, the component a*y greatly varies depending on a combination of j and i. Hence, in the case 
of estimating the noise X*, while changing Rj, only a plurality of relatively large values P in the component a*j are 
selected and only those i'-th ones of the components of the power P x . of Eq. 53 which correspond to the i-th compo- 
nents are selectively added to obtain P' Qp which is used as an evaluation measure-this increases the accuracy of esti- 

io mation of the value X y 

Furthermore, the estimation accuracy could be increased, for example, by a method in which a plurality of kinds of 
waveforms g(t), such as mentioned above, are prepared and evaluation values Px based on them are combined to 
adopt, for example, the minimum (or maximum) noise estimated value X*j. 

Although for the input signal and the signals taken in the process of noise detection, DFT and IDFT processing 
is have been used, instead r(t)xmsc{t) can be implemented by a convolution of Rj of Eq. 19 and C'j of Eq. 1 1 . Thus, the 
entire arithmetic processing in Figs. 1 , 3 and 4 can also be performed in the frequency domain. As a result, the matched 
filter output e l is obtained, which needs only to be compared with a preset threshold value. 

As described above, according to the present invention, the noise and signal components contained in the received 
frame signal can be detected and estimated by a simple method which multiplies the spread spectrum signal by a 
20 spreading signal twice, analyzes the re-spread output and solves simultaneous linear equations with L unknown varia- 
bles. 

Further, with a method which repeats an operation in which after the received signal is multiplied by an analyzing 
wave form, an excessively large or small amplitude component in the received signal spectrum is decreased or 
increased in such a manner as to reduce the power of that one of the frequency components of the multiplied output 

25 which intersects orthogonally a signal corresponding frequency component except noise, the noise component con- 
tained in the received frame signal can be estimated from the corrected value. 

By subtracting the noise component from the received frame signal, it is possible to obtain a received signal with 
no noise component or having the noise component reduced, so that information of the transmitted signal can be iden- 
tified by subjecting the received signal to ordinary demodulation processing. Except when the frequency characteristics 

30 of noise (the magnitude and phase polarity of the frequency spectrum) are closely resemble the frequency character- 
istics of the spreading code, noise can be estimated with high accuracy. Thus, the present invention enhances the SN 
ratio of the received signal and hence dramatically improves the error rate in the spread spectrum communication. 

Claims 

35 

1 . A method for separating and detecting a noise component contained in a spread spectrum signal that is performed 
at a receiver of a spread spectrum communication system is characterized in that: 

a spreading code sequences are used for corresponding to binary or multi-valued information, 
a received frame signal r(t), composed of a spreading code sequence m(t) of a desired station and noise n(t), 
is multiplied by an analyzing code sequence g(t) to obtain the multiplied output u(t); 

a set of N complex frequency components [UJ (where i=0, 1,2 N-1) obtained by discrete fourier transform 

analysis of the multiplied output u(t), 

the received frame signal with n(t)=0 is multiplied by g(t) to obtain the multiplied output y(t). 
a set of N complex frequency components pr] is obtained by discrete transform analysis of y(t); 
the set [U J of the complex frequency components is divided into the in-phase component [U,J and the quadra- 
ture phase component [U Qi ] to the set prj; 

the received frame signal r(t) is further analyzed to obtain a set [RJ of complex frequency components; 
the amplitude of a component R|< in the set of the complex frequency component [RJ, which has an excessively 
large or small amplitude as compared with a preset correcting spectrum curve, is corrected within a range in 
which a value, obtained by selectively adding together the powers of the set [UJ and the component [U Q J 
(where i=0.2 N-1) with respect to i, decreases or remains unchanged; 

the amplitude of a component R k , of a second excessively large or small amplitude is corrected within a range 
in which the power sum decreases or remains unchanged; 

these operations are repeated to reduce the power sum toward zero, then noise components X* and X' k are 
estimated from corrected values of the components l\ and R' k , and the estimated noise components are 
removed from the received frame signal r(t). 

2. The method of claim 1 for separating and detecting a noise component contained in a spread spectrum signal in a 
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spread spectrum communication system, characterized in that, of the frequency components of the received frame 
signal, only component projected to the phase of each frequency component [Cj obtained by a discrete Fourier 
transform analysis of the spreading code sequence m(t) of the desired station is extracted, and that the time signal 
obtained by the inverse discrete Fourier transform analysis is used as the received frame signal r(t). 

5 

3. The method of claim 1 for separating and detecting a noise component contained in a spread spectrum signal in a 
spread spectrum communication system, characterized in that a square wave spreading code sequence m(t) of the 
desired station or a sequence m L (t) obtained by band limiting it is used as an analyzing code sequence g(t). 

10 4. The method of claim 1 for separating and detecting a noise component contained in a spread spectrum signal in a 
spread spectrum communication system, characterized in that the analyzing code sequence g(t) is selected so that 
the multiplied output of the analyzing code sequence g(t) and the spreading code sequence m(t) of the desired sta- 
tion becomes a square wave that has a period kTc (where k=.... 1/3,1/2,1 ,2,3....). 

is 5. The method of claim 1 for separating and detecting a noise component contained in a spread spectrum signal in a 
spread spectrum communication system, characterized in that the analyzing code sequence g(t) is selected so that 
the multiplied output of the analysis code sequence g(t) and a sequence m L (t) obtained by band limiting the square 
wave spreading code sequence m(t) of the desired station is composed of a sine wave of a period kTc and a limited 
number of its high-frequency waves. 

20 

6. The method of claim 1 for separating and detecting a noise component contained in a spread spectrum signal, 
characterized by the use of a received signal containing noise and a time waveform obtained by limiting a spreading 
code sequence for use at the receiving side to a passband frequency band lower than a predetermined upper limit 
frequency. 

25 

7. The method to claims 1 through 5 for separating and detecting a noise component contained in a spread spectrum 
signal, characterized in that the multiplication of the received frame signal r(t) or multiplied output u(t) and the 
receiving side spreading code sequence g(t) in the time domain is performed by a convolution in the frequency 
domain. 

30 

8. The method of claims 1 through 5 for separating and detecting a noise component contained in a spread spectrum 
signal, characterized in that a transmitted signal is a time waveform obtained by the convolution of an impulse train 
produced by q-point sampling per chip time width of a square wave M-sequence and a sampling function waveform 
having its pass band limited to q/2, and that the receiver uses, as the despreading signal, the same waveform as 

35 the time waveform. 

9. An apparatus for separating and estimating a noise component contained in a spread spectrum signal in a receiver 
demodulating system in a spread spectrum communication system in which a spreading code sequence is made 
to correspond to binary or multi-valued information, 

40 CHARACTERIZED BY: 

means for despreading a received frame signal r(t) containing noise and the same spreading code sequence 
m(t) as that of a desired station; 

means for multiplying again the despread AC output, obtained by removing the DC component from the 
despread output u(t), with the spreading code sequence m(t) to obtain the re-spread output w(t); 

means for obtaining N respread frequency components Wi with i-times (where i=0,1,2 N-1) the frame rate 

that is the output obtained by a discrete Fourier transform analysis of the respread output w(t); 
means for solving simultaneous linear equation with N unknown variables consisting of a set of j-th (where 
j=0,1,2,..., N-1) noise frequency components Xj obtained by a discrete Fourier transform analysis of the 
received frame signal with the noise alone, and the coefficient matrix which relates the noise frequency com- 
ponents Xj and the respread frequency components Wi obtained in advance, and the set of Wj; and 
means for detecting the signal component in the frame obtained by removing the estimated noise wave form, 
inverse discrete fourier transform of the set Xj, from the received frame signal r(t). 
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